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In rotating machinery dynamics an orbit (Lissajous curve) represents the dynamic 
path of the shaft centerline motion during shaft rotation and resulting precession. 
The orbit can be observed with an oscilloscope connected to XY proximity probes. 

The orbits can also be simulated by a computer. 


OBJECTIVE 


The software for HP computer simulates orbits for two cases: 

(1) Symmetric orbit with four frequency components with different radial amplitudes 
and relative phase angles. 

(2) Nonsymmetric orbit with two frequency components with two different vertical/ 
horizontal amplitudes and two different relative phase angles. 

Each orbit carries a Keyphasor mark (one-per-turn reference). The frequencies, 
amplitudes, and phase angles, as well as number of time steps for orbit computa- 
tion, have to be chosen and introduced to the computer by the user. 

The orbit graphs can be observed on the computer screen. Hard copies are available 
from the corresponding plotter. 


RESULTS 

Figures 1 through 5 give the examples of the computer-generated orbits, simulating 
various rotating machine malfunctions. 
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QURDRflTURE 



-120E-0i -680e-02 +e06£+ee ^80E-0S +12BE-01 

DIRECT 


»**« DATA **** 


1. ) axx 

2. ) AYI 

3. ) AX2 
4.5 AY2 
3.> ALPX1 
€.5 ALPYi 
?.> ALPX2 

8. ) ALPY2 

9. > 0ME1 
18. > 0ME2 
11. > NMAX 


HORIZONTAL AMPLITUDE OF THE 1ST COMPONENT *************** 2 

VERTICAL AMPLITUDE OF THE 1ST COMPONENT ***************** 2 

HORIZONTAL AMPLITUDE OF THE 2ND COMPONENT *************** 9 

VERTICAL AMPLITUDE OF THE 2ND COMPONENT ***************** 9 

HORIZONTAL PHASE ANCLE OF THE 1ST COMPONENT CDEG3 ******* 98 

VERTICAL PHASE ANCLE OF THE 1ST COMPONENT CBEG3 ********* 80 

HORIZONTAL PHASE ANCLE OF THE 2ND COMPONENT CDEG3 ******* 0 

VERTICAL PHASE ANCLE OF THE 2ND COMPONENT CDEG3 ********* 8 

FREQUENCY OF THE 1ST COMPONENT CRAB^SECl **************** 188 

FREQUENCY OF- THE 2ND COMPONENT CRADsSECl **************** 43 

NUMBER OF 'TIME '-STEPS ********************************** 933 


»*** ROTOR ORBITAL MOTION **** 


x» < Rx O c os < o**g*i ft phax 1 > + < Ax2> cos < om*q*Z*z +*1 phax2 > 
y*< Ay 1 >si n < on»*gal *% ph ay 1 > ♦ < Ay 2 ) s i n < o«*g *2*t ph ay 2 ) 


Figure 1, - Orbit simulating oil whirl. 
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QUHDRRTURE 



-ia*£-*i -i«ee*S2 *aeac*e» «*mc m ♦i2«c-n 

DIRECT 


**** DATA **** 


L) ftXl 

2 . ) ayi 

3. > AX2 

4. ) AY2 
3.) ALPX1 

5. > ALPY1 

7.3 RLPX2 
S.> ALPY2 

3.3 0M61 

19.3 0ME2 

11. 3 NMAX 


HORIZONTAL AMPLITUDE OF THE 1ST COMPONENT *************** 2 

VERTICAL AMPLITUDE OF THE 1ST COMPONENT ***************** 2 

HORIZONTAL AMPLITUDE OF THE 2ND COMPONENT *************** 3 

VERTICAL AMPLITUDE OF THE 2ND COMPONENT ***************** 9 

HORIZONTAL PHASE ANGLE OF THE 1ST COMPONENT CBEG3 ******** 223 
VERTICAL PHASE ANGLE OF THE 1ST COMPONENT CBEC3 ********* 223 

HORIZONTAL PHASE ANGLE OF THE 2ND COMPONENT CDEG3 ******* 20 

VERTICAL PHASE ANGLE OF THE 2ND COMPONENT C3EG3 ********* 200 

FREQUENCY QF THE 1ST COMPONENT CRAD^S EC3 **************** 100 
FREQUENCY OF THE 2ND COMPONENT CRAD/SEC3 **************** *30 
NUMBER OF 'TIME'- STEPS ********************************** 300 


**** ROTOR ORBITAL MOTION **** 


x*< Ax 1 > c os < o»*gsl *a ! ph<axl > +< Ax23c os < oaicg phax2 ) 

y»< Ayl3si n<o«*gal+*!phayl 3+<Ay23*in<oift«ga2+*! phay23 

Figure 2. - Orbit simulating vertical partial rotor-to-stator rub. 



QUADRATURE 



-laee-ei -600E-02 +000^^0 ^seae-ea +iaaE-ai 

DIRECT 


**** DATA **** 


1. ) AX1 

2. ) AY1 
3*) AX2 
4,) AY2 

3. ) ALPXl 
Sc) ALPY1 
7c > ALPX2 
8.) ALPY2 
9») 0ME1 

18.) 0ME2 
11.) NMAX 


HORIZONTAL AMPLITUDE OF THE 1ST COMPONENT *************** 2 

VERTICAL AMPLITUDE OF THE 1ST COMPONENT ***************** 2 

HORIZONTAL AMPLITUDE OF THE 2ND COMPONENT *************** 9 

VERTICAL AMPLITUDE OF THE 2ND COMPONENT ***************** 1 

HORIZONTAL PHASE ANGLE OF THE 1ST COMPONENT CDEG3 ******* 0 

VERTICAL PHASE ANGLE OF THE 1ST COMPONENT EDEG3 ********* 0 

HORIZONTAL PHASE ANGLE OF THE 2ND COMPONENT CDEC3 ******* 0 

VERTICAL PHASE ANGLE OF THE 2ND COMPONENT EDEG3 ********* 90 

FREQUENCY OF THE 1ST COMPONENT CRAD^SEC3 **************** 100 
FREQUENCY OF THE 2ND COMPONENT CRAD^SEC3 **************** -30 
NUMBER OF * TIME' -STEPS ********************************** 300 


**** ROTOR ORBITAL MOTION **** 


x« < Ax 1 >c os < ora *g*l *v+&1 phax 1 ) +< Ax2 ) c os < o»«ga2*t +*1 phax2 ) 
y*<Ayl )si n<o»*gal*i+alphayl )+<Ay2)si n<on*ga2*t+al phay2) 

Figure 3 . - Orbit simulating horizontal partial rotor-to-stator rub. 
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QUADRATURE 



DIRECT 


x 


■ ClB3EXPejCS0t+7033 «■ C183EXPC JC1B0X+B33 + C 13EXPC JC 130«*a3 3 
| +• a)EXP{j<20et-H3)) 

FIRST COMPONENT IS SYNCHRONOUS 


ROTOR ORBITRL MOTION 


Figure 4. - Orbit simulating cracked shaft. 
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QUADRATURE 


+250E-02 


+125E-02 


+000E+00 


12SE-02 



DIRECT 


**** BATA **** 


1. ) axi 

2. > AYl 
3* ) AX2 
4c > AY2 
3c > ALPX1 
So) ALPY1 
?.) ALPX2 

3. > ALPY2 
9* > 0ME1 

13c) 0ME2 
11.) NMAX 


HORIZONTAL AMPLITUDE OF THE 1ST COMPONENT *************** -2 
VERTICAL AMPLITUDE OF THE 1ST COMPONENT ***************** - e 23 
HORIZONTAL AMPLITUDE OF THE 2ND COMPONENT *************** .5 

VERTICAL AMPLITUDE OF THE 2ND COMPONENT ***************** i 
HORIZONTAL PHASE ANGLE OF THE 1ST COMPONENT CDEGI ******* 3 

VERTICAL PHASE ANGLE OF THE 1ST COMPONENT C DEG 3 ********* 133 

HORIZONTAL PHASE ANGLE OF THE 2ND COMPONENT CBEG3 ******* 3 

VERTICAL PHASE ANGLE OF THE 2ND COMPONENT CDEG3 ********* 133 

FREQUENCY OF THE 1ST COMPONENT CRAD^SECI **************** 133 

FREQUENCY OF THE 2ND COMPONENT CRAD^SECI **************** 97 

NUMBER OF . 'TIME' -STEPS ********************************** 933 


•*** ROTOR ORBITAL MOTION **** 


x* < Ax 1 ) c os < om*g &1 *t +* 1 ph 1 ) * < Ax2 ) cos < o»«g*2* t +*1 phax2 > 
y»<Ayi)sin< ora e g *1 *t ** ? ph ay 1 ) + < Ay 2 ) s i n < o»*g *2 * t ** I phay 2 ) 


Figure 5. - Orbit simulating rotor instability. 
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